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Abstract 

Chaos degree defined through two complexities in information dy- 
namics is applied to some deterministic dynamical models. It is shown 
that this degree well describes the chaostic feature of the models. 

1 Introduction 

There exist several approaches in the study of chaotic behavior of dynam- 
ical systems using the concepts such as entropy, complexity, chaos, fractal, 
stochasticity [1], [| £|, |5], |], [U|. In 1991, one of the authors proposed Infor- 
mation Dynamics (ID for short) [[HJ [TlJ to treat such chaotic behavior of 



systems altogether. ID is a synthesis of the dynamics of state change and 
the complexity of systems, and it is applied to several different fields such as 
quantum physics, fractal theory, quantum information and genetics^. 

A quantity measuring chaos in dynamical systems was defined by means 
of two complexities in ID, and it is called chaos degree. In particular, among 
several chaos degrees, an entropic chaos degree was introduced in [0, ||, and 
is applied to logistic map to study its chaotic behavior. This chaos degree 
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has several merits compared with usual measures of chaos such as Lyapunov 
exponent. 

In Section2, we review the complexity in ID and the chaos degree (CD 
for short). In Section 3, we remind the entropic chaos degree and the Lya- 
punov exponent (LE for short). In Section 4, the algorithm computing the 
entropic chaos degree is shown. In Section 5, we compute CD and LE for 
Bernoulli shift, Baker's transformation and Tinkerbell's map, then we discuss 
the merits of the entropic chaos degree. 

2 Complexity of Information Dynamics and 
Chaos Degree 

Information dynamics provides a frame to study the state change and the 
complexity associated with a dynamical system. We briefly explain the con- 
cept of the complexity of ID in a bit simplified version (see||ll||). 

Let (A,&,a(G)) be an input (or initial) system and (A, &,a(G)) be an 
output (or final) system. Here A is the set of all objects to be observed and 
& is the set of all means for measurement of A, a(G) is a certain evolution 
of system. Often we have A = A, (5 = 6, a = a. 

For instance, when A is the set M(Q) of all measurable functions on a 
measurable space (Q, T) and &(A) is the set -P(fi) of all probability measures 
on Q, we have usual probability theory, by which the classical dynamical 
system is described. When A is the set B{H) of all bounded linear operators 
on a Hilbert space TC and &(A) is the set &(H) of all density operators on 
7i, we have a quantum dynamical system. 

Once an input and an output systems are set, the situation of the input 
system is described by a state, an element of & , and the change of the state 
is expressed by a mapping from & to &, called a channel. The concept of 
channel is fundamental both in physics and mathematics 0. Moreover, there 
exist two complexities in ID, which are axiomatically given as follows: 

Let (At, ©t, a*(Cr*)) be the total system of (A,&,a) and (A,&,a), and 
let C (if) be the complexity of a state if and T (if] A*) be the transmitted 
complexity associated with the state change if — > A*ip. These complexities C 
and T are the quantities satisfying the following conditions: 

(i) For any f G (5, 

C(<f) > 0, T(ip; A*) > 0. 
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(ii) For any orthogonal bijection j : ex&^ex& ( the set of all extreme 
points in & ), 

T(j(f);A*)=T(if;A*). 

(hi) For $ = Lp g) -0 e ©t, 

(7(*) = C(^) + C7(Vi). 

(iv) For any state and a channel A*, 

< T{tp\ A*) < (%). 

(v) For the identity map "id" from & to &. 

T(<p;id) = C( ( p). 

When a state if changes to the state A* if, a chaos degree (CD) [|11|] w.r.t. 
if and A* is given by 

D(if;A*) = C(A*f)-T(if;A*). 

Using the above CD, we observe chaos of a dynamical system as 

CD > <^=^ chaotic. 
CD = stable. 

3 Entropic Chaos Degree and Lyapunov Ex- 
ponent 

Chaos degree in ID was applied to a smooth map on R and it is shown 
that this degree enables to describe the chaotic aspects of a logistic map as 
well|| [12]. Here we briefly review the chaos degree defined through classical 
entropies. For an input state described by a probability distribution p = (pj) 
and the joint distribution r = (rij) between p and the output state p = 
A*p = (pj through a channel A*, the Shannon entropy S (p) — — J2iPi^°SPi 
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and the mutual entropy I(p;A*) = J2ij r i,j^°S'^p' satisfy all conditions of 
the complexities C and T, then the entropic chaos degree is defined by 



D(p;A*) = C(A*p)-T(p;A*) 
= S(p)-I(p;A*), 

This entropic chaos degree is nothing but the conditional entropy of apos- 
teriori state w.r.t. the channel A*. The characteristic point of the entropic 
chaos degree is easy to get the probability distribution of the orbit for a 
deterministic dynamics, which is discussed in the next section. 

Lyapunov exponent (LE) is used to study chaotic behavior of a determin- 
istic dynamics. The Lyapunov exponent A (/) for a smooth map / on R is 
defined by [0: 



n 



fc=i 



— (x (k) 
dx 



A(/)= lim A n (/) A n (/) 

where x^ 1 ' = f (a^™" 1 )) for any n G N. 

For a smooth map / = (/i, • • - , f m ) on R m , the vector version of LE is 
defined as follows: Let x° be an initial point of R m and x^ = f (x^ ri_1 - ) ) 
for any n G N. After n times iterations of / to x ' the Jacobi matrix J n (x°) 



of x {n) 



x 



(n) 



(n) 



w.r.t. x^is 
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Then the Lyapunov exponent A (/) of x^ is defined by 
AC/ - ) = log/ii,/Zjfc = lim (/x^)" (fc= 1,2,- • 



,m . 



Here ^x)! is the fcth largest square root of the eigenvalues of the matrix 

J n (*<°>) J n (x(°)) T 

An orbit of the dynamical system described by / is said to be chaotic 
when the exponent A (/) is positive, and to be stable when the exponent is 
negative. The positive exponent means that the orbit is very sensitive to the 
initial value, so that it describes a chaotic behavior. Lyapunov exponent is 
difficult to compute for some models (e.g., Tinkerbell map) and its negative 
value is not clearly explained. 
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4 Algorithm for Computation of Entropic 
Chaos Degree 

It is proved || that if a piecewise monotone mapping / from [a, b] m to [a, b] m 
has non-positive Schwarzian derivatives and does not have a stable and pe- 
riodic orbit, then there exists an ergodic measure /i on the Borel set & of 
[a, b] m , absolutely continuous w.r.t. the Lebesgue measure. 
Take a finite partition {Ak} of I = [a, b] m such as 

I = \jA k (^f|^ = 0,*^j). 

k 

Let 15*1 be the number of the elements in a set S. Suppose that n is sufficiently 

large natural number and m is a fixed natural number. Let p^ = [P^j be 

the probability distribution of the orbit up to n -th step, that is, how many 
x (k) ^ _ m _|_ . . . ^ m _|_ n ^ are m j^.. 

(n) _ | G N; x (k) E Ai, m < k < m + n}\ 



It is shown that the n — ► oo limit of p\ n exists and equals to n {Ai). The 
channel A* is a map given by p( n+1 ) = A*p ( - n \ Further, the joint distribution 
r (n,n+i) _ ^ r ^. ,n+1 ^f r a sufficient large n is approximated as 

(n,n+i) _ |{fceN; (x^, f (x^)) e Ai x Aj , m < k < m + n}\ 
T%J = n ' 

Then the entropic chaos degree is computed as 



£>(p(");A*) = S (p in+1) ) - I (p (n) ;A*) 

= -ep. ( " +i, ^p. ( " + " - s>-r +i, '° g ^y 

5 Entropic Chaos Degree for Some Determi- 
nistic Dynamical Models 

In this section, we study the chaotic behavior of several well-known deter- 
ministic maps by the entropic chaos degree. 
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(1) Bernoulli shift 

Let / a map from [0, 1] to itself such that 



(5.1) 



2ax^ (0 < iW < 0.5) 

a (2x^ - 1) (0.5 < x^ < 1) ' 

where x^ G [0,1] and < o < 1. 

Let us compute the Lyapunov exponent and the entropic chaos degree 
(ECD for short) for the above Bernoulli shift / . 

The orbit of the equation (5.1) is shown in Fig. 1: 




0.1 0.2 0.3 0.1 0.5 0.6 0.7 0.8 0.9 1 

Figl. The bifurcation diagram for Bernoulli shift 

The Lyapunov exponent A n (/) is log 2a for the Bernoulli shift (Fig. 2). 




Fig. 2 Lyapunov exponent of Bernoulli shift 

On the other hand, the entropic chaos degree of the Bernoulli shift is 
shown in Fig. 3. 
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Here we took 740 different a's between and 1 with 



n 



i i + l 
2000' 2000 _ 
= 100000. 



(),•••, 1999) 



(2) Baker's transformation 

We apply the chaos degree to a smooth map on R 2 . Let us compute the 
Lyapunov exponent and the ECD for the following Baker's transformation 

fa ■ 

fa ^) = fa (^ x l \ x 2 ^ 

2ax[ n) ia4 n) ) (o < xf ] < 0.5^ 




[a (2xf ] - l) , \a (4 n) + l)) (0.5 < x^ < 1 
where (x^ , x^) E [0, 1] x [0, 1] and < a < 1 
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The orbit for each a is shown in Fig. 4. ~9. 




0.2 OA 0.6 o.; 

Fig.4 Orbit of /. w.r.t. o < a < 0.5 




02 OA 0.6 0.8 

Fig.5 Orbit of /„ w.r.t. a = 0.6 




02 OA 0.6 Of. 

Fig.6 Orbit of /„ w.r.t. a = 0.7 
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02 OA 0.6 0.8 

Fig.7 Orbit of/, w.r.t. a = 0.8 
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Fig.8 Orbit of /„ w.r.t. a = 0.9 




Fig.9 Orbit of /„ w.r.t. a = 1.0 



These figures show that the larger a is, the more complicated the orbit is. 
The maximum Lyapunov exponent A* (/) is log 2a for Baker's transformation 
(Fig.10). 
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Mf) 




Fig. 10 Lyapunov exponent of Baker's transformation 

On the other hand, the ECD of Baker's transformation is shown in Fig. 11. 
ECD 




OjO 0.1 0:2 0.3 0.4 0J5 0j6 0.7 OjS 0.-3 ijO 

Fig. 11 ECD of Baker's transformation 



Here we took 740 different a's between and 1 with 



A 



i i + 1 

Too' ioo 



x 



100' 100 



(U = 0,--.,99) 



n = 100000. 



(3) Tinkerbell map 



Let us compute the CD for the following two type Tinkerbell maps f a and 
f b on 1= [-1.2,0.4] x [-0.7,0.3]. 
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fa(x^) = / a (4 B U B) ) 

= ((^! n) ) 2 - (4 n) ) 2 + ^ + c 2 xi n \2x^x^ + c 3 x[ n) + c 4 x { ^ , 
h(x {n) ) = f b (x[ n \xP) 

= ((^S n) ) 2 - (^2 n) ) 2 + cixS n) + c 2 4 n) , 2x^x { ^ + bx™ + c 4 x { ^ , 

where (x^^^ e 7,-0.4 < a < 0.9,1.9 < b < 2.9, (c u c 2 , c 3 , c 4 ) = 

(-0.3, -0.6, 2.0, 0.5) , and (x^ , 4 0) ) = (0.1, 0.1) . 

Let us plot the points (x^, ) f° r 3000 different n's between 1001 and 
4000. 

In stable domain, the number of the points ^x[ n \ x^ j is finite because 

the point (x^jX^ 1 ^ periodically appears in time n. Fig. 12 and Fig. 13 are 
examples of the orbits of f a and in a stable domain. 



0.6 
OA 
' 0.2 
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-1JD — J8 -0.6 -0.4 -02 Oi * 0.2 0.1 
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Fig. 12 Orbit off,, w.r.t. a =0.243 
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Fig. 1 3 Orbit of /„ w.r.t. b = 2.65 



On the other hand, the point {x^^x^^j take random value in time n in 
chaotic domain. Fig.514 and Fig. 15 are examples of the orbits of f a and /& 
in a chaotic domain. 




Fig. 14 Orbit of/, w.r.t. a =0.670 




Fig. 15 Orbit of f„ w.r.t. b = 2.8 



The ECD of Tinkerbell map f a and f b are shown in Fig. 16 and Fig. 17. 
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ECD 




-02 -0.1 0.1 0.2^0.3 0.4 0.5 0.6 0.7 0.8 

Figia 

Fig. 16 ECD for Tinkerbell map/, 
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Fig. 15 

1 





2.0 2.1 2.2 2.3 2.4 2.5 2.6^2.7 2.8 ' 2.9 3.0 

Fig 13 

Fig. 17 ECD for Tinkerbell map/, 



a 



Here we took 740 different a's between —1.2 and 0.9 and 740 different 6's 
between 1.9 and 2.9 with 



% i + 1 
100' 100 



x 



3 3 + 1 



n = 



100' 100 

(i = -120, -119, • • • , -1, 0, 1, • • • , 38, 39) 
(j = -70, -69,-.., -1,0,1,.-., 28, 29) 
100000. 
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6 Conclusion 



From our results, when the orbits for Bernoulli shift, Baker's transformation 
or Tinkerbell map are chaotic, both Lyapunov exponent and chaos degree are 
positive. However, our chaos degree can resolve some inconvenient properties 
of the Lyapunov exponent in the following senses: 

(1) Lyapunov exponent takes negative and sometimes — oo. For instance, 
although the exponent of Bernoulli shift and Baker's transformation 
can not be defined for a = 0, the ECD is always positive and defined 
for any a > 0. 

(2) It is difficult to compute the Lyapunov exponent for the Tinkerbell 
maps f a and ft, because it is difficult to compute /" and for large 
n. On the other hand, the ECD of f a and f are easily computed. 

(3) Generally, the algorithm for CD is much easier than that for Lyapunov 
exponent. 
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